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Abstract 

A monomial dynamical system / : over a finite field i^T is a nonlinear deterministic time discrete 

dynamical system with the property that each component function fi : K"' ^ if is a monic nonzero monomial 
function. In this paper we provide an algebraic and graph theoretic framework to study the dynamic 
properties of monomial dynamical systems over a finite field. Within this framework, characterization 
theorems for fixed point systems (systems in which all trajectories end in steady states) are proved. In 
particular, we present an algorithm of polynomial complexity to test whether a given monomial dynamical 
system over a finite field is a fixed point system. Furthermore, theorems that complement previous work are 
presented and alternative proofs to previous results are supplied. 



1 Introduction 

Time discrete dynamical systems over a finite set X are an important subject of active mathematical 
research. One relevant example of such systems are cellular automata, first introduced in the late 
1940s by John von Neumann (e.g., [2|). More general examples of time discrete dynamical systems 
over a finite set X are non- deterministic finite state automata (e.g., [20]) and sequential dynamical 
systems [3]. 

Deterministic time discrete dynamical systems over a finite field are mappings / : K"" — >■ K"", 
where if is a finite field and n € N the dimension of the system. They constitute a particular class 
of deterministic time discrete dynamical systems over a finite set X, namely, the class in which the 
finite set X can be endowed with the algebraic structure of a finite field. This property allows for 
a richer mathematical framework within which these systems can be studied. For instance, it can 
be shown that every component function fi : K"^ — >■ if is a polynomial function of bounded degree 
in n variables (see, for example, pages 368-369 in [16] or 3.1 in [8j). 

The study of dynamical systems generally addresses the question of the system's long term 
behavior, in particular, the existence of fixed points and (limit) cyclic trajectories. (The state of 
the system evolves by iteration of the function / starting from given initial conditions xq G if".) In 
this paper we provide an algebraic and graph theoretic framework to study a very specific class of 
nonlinear time discrete dynamical systems over a finite field, namely, monomial dynamical systems 
over a finite field. In such systems, every component function fi : if" ^ if is a monic nonzero 
monomial function. 
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Some types of monomial systems and their dynamic behavior have been studied before: mono- 
mial cellular automata [H], [1], Boolean monomial systems [6], monomial systems over the p-adic 
numbers [15], [18] and monomial systems over a finite field [21], [1], [5]. [6| proved a necessary 
and sufficient condition for Boolean monomial systems to be fixed point systems (systems in which 
all trajectories end in steady statesj^. This condition could be algorithmically exploited. Indeed, 
the authors make some suggestive comments in that direction (see 4.3 in [6]). Moreover, the pa- 
per describes the structure of the limit cycles of a special type of Boolean monomial systems. [5] 
presents a necessary and sufficient condition for monomial systems over a finite field to be fixed 
point systems. However, this condition is not easily verifiable and therefore the theorem does not 
yield a tractable algorithm in a straightforward way. 

Our work was strongly influenced by [6], [5] and [1]. However, we took a slightly different 
approach. The mathematical formalism we developed allows for a deeper understanding of mono- 
mial dynamical systems over a finite field. In particular, we present an algorithm of polynomial 
complexity to test whether a given monomial dynamical system over a finite field is a fixed point 
system. Furthermore, we obtain additional theorems that complement the work of [6], [^ and 
provide alternative proofs to many results in [6]. Our formalism also constitutes a basis for the 
study of monomial control systems, to be presented elsewhere. 

It is pertinent to mention the work of [9] regarding linear time discrete dynamical systems over a 
finite field, in which the number of limit cycles and their lengths is linked to the factorization (in so 
called elementary divisor polynomials) of the characteristic polynomial of the matrix representing 
the system. (See also [TIj for a more mathematical exposition and [20], [19] for applications of the 
Boolean case in control theory.) Furthermore, in [T7|, the affine case (a linear map followed by a 
translation) was studied. An interesting contribution was made by Paul Cull ([7]), who extended the 
considerations to nonlinear functions, and showed how to reduce them to the linear case. However, 
Cull's approach does not yield an algorithm of polynomial complexity to solve the steady state 
system problem. Moreover, according to [T2|, this might in general not be possible as a matter of 
principle. 

The organization of this article is the following: 

Section 2 establishes an algebraic and graph theoretic framework within which monomial dy- 
namical systems over a finite field are studied. It starts with some basic definitions and algebraic 
results (some of which are proved in the appendix) and leads the reader to the first important 
result: Theorem 2, which states that the monoid of n-dimensional monomial dynamical systems 
over a finite field K is isomorphic to a certain monoid of matrices. Section 2 finishes with propo- 
sitions about the relationship between the matrix F corresponding to a monomial system / (via 
the isomorphism mentioned above) and the adjacency matrix of the dependency graph of / (to be 
defined below). 

Section 3 is devoted to the characterization of fixed point systems. These characterizations are 
stated in terms of connectedness properties of the dependency graph. We provide some necessary 
and sufficient conditions for a system to be a fixed point system (Theorems 6 and 8). Moreover, we 
prove several sufficient conditions for special classes of monomial dynamical systems over a finite 
field K. 

Section 4 presents an algorithm of polynomial complexity to test whether a given monomial 
dynamical system over a finite field K is a fixed point system. A detailed complexity analysis of 
the algorithm is provided. 



^ This problem is refered to as the steady state system problem, see [12] , 
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2 Algebraic and graph theoretic formalism 

In this section we will introduce the monoid of n-dimensional monomial dynamical systems over 
a finite field F^. Furthermore we will show that this monoid is isomorphic to a certain monoid of 
matrices. This result establishes that the composition f o g oi two monomial dynamical systems 
/, g is completely captured by the product F ■ G oi their corresponding matrices. In addition, 
we will introduce the concept of dependency graph of a monomial dynamical system / and prove 
that the adjacency matrix of the dependency graph is precisely the matrix F associated with / 
via the isomorphism mentioned above. This finding allows us to link topological properties of the 
dependency graph with the dynamics of /. 

Definition 1 (Notational Definition) Since for every finite field K there is a prime number 
p G N f' the characteristic of K) and a natural number n G N such that for the number of elements 
\K\ of K it holds 

\K\ 

we will denote a finite field with ¥q, where q stands for the number of elements of the field. It is of 
course understood that q is a power of the (prime) characteristic of the field. 

Definition 2 Let Fg be a finite field. The set 

:={0,...(?-l} CNo 

is called the exponents set to the field Fg. 

Definition 3 Let Fg be a finite field. A map f : F" — > F" is called a monomial dynamical system 
over Fg if for every i E {1, ...,n} there exists a tuple {Fn, ...,Fin) G -Eg such that 

= ^..x^- V X G F^ 

Remark 4 As opposed to JBjl, we exclude in the definition of monomial dynamical system the 
possibility that one of the functions fi is equal to the zero function. However, in contrast to we 
do allow the case fi = l in our definition. This is not a loss of generality because of the following: 
If we were studying a dynamical system f : Fg — > Fg where one of the functions, say fj, was equal 
to zero then for every initial state x € F^ after one iteration the system would be in a state f{x) 
whose jth entry is zero. In all subsequent iterations the value of the jth entry would remain zero. 
As a consequence, the long term dynamics of the system are reflected in the projection 

T^jiy) ■= {yii--,yj-i,yj+i,-;ynY 

and it is sufficient to study the system 

I fi{yi, ■■■,yj~i,0,yj+i, ...,yn) \ 

^ fj-i{yi,--,yj-i,0,yj+i,...,yn) 
fj+i{yi, ■■■,yj~i,o,yj+i, ■■■,yn) 

\ fniyi, ■■■,yj-i,0,yj+i, ...,yn) ) 
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In general, this system f could contain component functions equal to the zero function, since every 
component fi that depends on the variable xj would become zero. As a consequence, the procedure 
described above needs to be applied several times until the lower n' -dimensional system obtained 
does not contain component functions equal to zero. The long term dynamics of f are reflected in 
the projection to an n' -dimensional subspace, in particular, all the cycles and fixed points of f are 
located in this lower dimensional space. Moreover, points located outside this lower dimensional 
subspace are transient states of the system. It is also possible that this repeated procedure yields 
the one dimensional zero function. In this case, we can conclude that the original system f is a 
fixed point system with (0, 0) G as its unique fixed point. Note that this procedure reduces the 
dimension by s where < s < n. As a consequence, the procedure needs to be iterated at most n 
times. 

As stated in Theorem 16 and Theorem 20 of [S], every function h : — > Fg is a polynomial 
function in n variables where no variable appears to a power higher or equal to q. Calculating the 
composition of a dynamical system / : F^ F^ with itself, we face the situation where some of 
the exponents exceed the value q — 1 and need to be reduced according to the well-known rule 



This process can be accomplished systematically if we look at the power x^ (where p > q) as a 
polynomial in the ring Fg[r] as described in the Lemma and Definition below. But first we need 
an auxiliary result: 

Lemma 5 Let Fq be a finite field and a gNq a nonnegative integer. Then 



Proof. If a = X{q - 1) then x" = jj^C^-i) = (x^^^^))^ = 1 V x € Fg\{0} by ([I]). Now assume 
j;" = 1 V X G Fg\{0} and write a = a{q — 1) + s with suitable a € No and < s < (g — 1). Then it 
follows 

1 = x" = = rE^('?-i)x^ = x'\/ xe Fg\{0} 

As a consequence, the polynomial — r" € Fg[r] has 

\Fq\ -l = q-l>s = deg(r' - r) 

roots in Fq and must be therefore of degree s = q — 1. Thus a = (a + l){q — 1). ■ 

Lemma 6 (and Definition) Let Fq be a finite field and c G No a nonnegative integer. The degree 
of the (unique) remainder of the polynomial division r'^ ^ (r"? — r) is called redq{c). redq{c) satisfies 
the following properties 

1. redq{redq{c)) = redq{c) 

2. redq{c) = <^ c = 

3. For a, 6 G No, x'' = V X G Fg <^ redq{a) = redq(b) 



= a V a G Fg 



(1) 



x'' = 1 V X G Fq\{0} 4^ 3 A G No : a 



\{q-l) 



4. For a, 6 G N, redq{a) = redq{b) -^Ba&Z: a = b + a{q — 1) 
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Proof. By the division algorithm there are unique g,r ^ ^gl''"] with either r = or deg(r) < 
deg(r^ — r) such that 

= gij'^ — t) +r 

If we look at the corresponding polynomial function^ defined on Fg it follows by ([1]) 

x'^ = r{x) V X e Fg (2) 

In particular, r 7^ 0. From the division process it is also clear that r must be a monomial and we 
conclude r = r^'^'^i^^^ with redq{c) < q. The first property follows trivially from the fact redq{c) < q. 
The second property follows immediately from evaluating the equation x" = (i.e. equation 

([2])) at the value x = 0. The third property is shown as follows: By the division algorithm 3i ga 
, gb , ra , n £ Fg[r] such that 

= ffa(r''-r)+r, = ff,(r''-r)+r-'^'(") (3) 
r' = 56(r''-T)+r, = 56(r«-r)+r-'^'('') 

From = V X € Fg now we have 

^red,{a) ^ ^red,(b) ^ 

and since redq{a),redq{b) < q we get redq{a) = redq{b). On the other hand, from redq{a) = redq{b) 
it would follow from equations ([3]) 

- gair'^ - r) = t'^ - gbir'^ - t) 

and thus by ([I]) 

x'^ = x^ V X € Fq 

Last we prove the fourth claim: If redq{a) = redq{b) then by 3. we have 

x'^ = x'' V X € Fq 

Now assume wlog a >b and d := a — 6 G Nq. Then the last equation can be written as 

x^x'^ = X* V X G Fg 

yielding 

x'^ = 1 V X G Fq\{0} 

By Lemma [5] we have 3 a G No : = a{q — 1) and therefore a = b + a{q — 1) or 6 = a — a{q — 1). 
Now assume the converse, namely 3aGZ:a = 6 + a{q — 1). Assume wlog a > (otherwise 
consider b = a — a{q — 1)). Then we would have 

T = T ^- 'T 



If r G Fq[r] is a polynomial of degree n, i.e. r = ^ liT*, then r is defined as the polynomial function 

i=0 



i=0 
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and thus by Lemma [5] 

rr'^ = X* V X G Fg\{0} 

Since a, 6 > we also have 



Remark 7 From the properties above we have = x^'^'^''^"'^ V x € Fg. 

The "exponents arithmetic" needed when calculating the composition of dynamical systems 
/, g : Fg Fg can be formalized based on the reduction algorithm described by the previous 
lemma. Indeed, the set 

Eg = {0,l,...,{q-2),{q-l)}cZ 

together with the operations of addition a © 6 := redq{a + h) and multiplication a»b := redq{ab) is 
a commutative semiring with identity 1. We call this commutative semiring the exponents semiring 
of the field Fg. This result is proved in the Appendix (see Theorem [62]) . We also defer to the 
appendix the proof of the following lemma: 

Lemma 8 Let n € N 6e a natural number, Fq be a finite field and Eq the exponents semiring of 
Fq. The set M{n x n; Eq) ofnxn quadratic matrices with entries in the semiring Eq together with 
the operation ■ of matrix multiplication ( which is defined in terms of the operations and • on the 
matrix entries) over Eq is a monoid. 

Remark 9 (and Definition) The operation redq ■.'No^Eq can be extended to matrices M{nxn; 
No) by applying redq to the entries of the matrix. We call this extension mredq : M{n x n; 
No) — > M(n x n; Eq). See Remark^U^ in the appendix for further details. One important property 
of mredq shown in Remark is 



mredq{A) = A = Q (4) 
Definition 10 Let Fg be a finite field and n, m E N natural numbers. The set 

MF:i,{¥q) : ={/:F^^F^| 

3F G M(n X m;Eg) : /i(x) := xfi-.x^^™ V X G F^} 

is called the set of n- dimensional monomial mappings in m variables. 

Lemma 11 Let Fg be a finite field and n,m,r £ N natural numbers. Furthermore, let f £ 
MF^{Fq) and g G MF^{Fq) with 

fi{x) = x^'\..Xn'" y X eFq, i = l,...,m 
gj{x) = xf^^-.x™^™ V X G F^, j = 1, ...,r 

where F G M{m x n; Eq) and G G M{r x m; Eq). Then for their composition g o f : F^ — > F^ it 
holds 

m 

igof)kix)=llx, '=1 VxGF-, ke{l,...,r} 

i=i 
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Proof. From the definition It follows for every k € {1, ...,r} 



{go f)k{x) = ll{Mx)f^' = ll{llx; 

1=1 1=1 j=i 

For a fixed but arbitrary m G N we will prove the claim using induction on the dimension n of 

_m_ 



m „ m 



the mapping g o f. For n = 1 we have {g o f)k{x) = = Yl xf'''^^'^ = x 

1=1 1=1 

m 

redgiY^GkiFn) 

Xi (see Remark [7]), thus the claim holds in dimension 1. Now we consider the case 

n + 1 : 



m n+1 

kl 



{9of)k{x) = \{C[{xf'f' 

1=1 j=i 



and by induction hypothesis 



n(<:'«)"n<")°" 

1=1 j=i 

ml n 

1=1 \ j=i 

m m n 

1=1 1=1 j=i 



E^Gfci-F'i(n+i) " redq{Y,GuFij) 
^(n+1) ii^J 

m m 

^(n+1) li^i'"' 



n+1 

n 

i=i 



7- • i = l 



Remark 12 (and Lemma) // we generalize the matrix multiplication defined on the monoid 
M{n X n; Eg) for matrices F G M(m x n; Eg) and G G M(n x m; £^q) i/ien we can write 

n 

{gof)k{x) = J{x,^''-^^>^^^x^Y^g, fcG{l,...,n} 

3=1 
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To see this, apply the Lemmas [61\ andlTl\ as well as the definitions of Q) and » to Yl Xj^^'^^^i 

n n 
3=1 3=1 



m ^mj ) 



"3 
3=1 
n 

— ^xedg{redg{GkiFij)+...+redg{GkmFmj)) 

m 

" redgiY^GkiFi^) 
= U^3 - 

i=i 

= {g°f)k{x) 

Theorem 13 Let Fq be a finite field. The set 

MF„"(F,) : = {/ : ^ I 

3F G M(n X : := xf"...x^^'" V X G F^} 

of all monomial dynamical systems overFq together with the composition o of mappings is a monoid. 



Proof. By Lemma [TT] the set MF^{Fq) is closed under composition. Composition of mappings 
is trivially associative. The identity function 

Id : F^ ^ F^ 

X ^ X 

is a monomial function and is therefore the identity element of the monoid {MF^(Fq), o). ■ 

Theorem 14 The monoids M{n x n; Eq) and MF^(Fq) are isomorphic. 
Proof. From the definition of MF^{Fq) it is clear that the mapping 

^ : M{n X n; Eg) MF^(Fg) 
G ^ -^{G) 



such that 



^'(G)i(x) := Xi'^...Xn'" for i = 1, ...,n 



is a bijection. Moreover, ^{I) = id. In addition, by Remark 1121 it follows easily that 

^{F-G) = ^'(F)o^(G) 



Remark 15 (and Definition) For a given monomial dynamical system f G MF^(Fq) the matrix 
F := '^~^{f) is called the corresponding matrix of the system f. For a matrix power in the monoid 
M{n X n; Eq) we use the notation F ™. By induction it can he easily shown 
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Remark 16 (and Definition) The image of the n x n zero matrix S M(n x n; Eq) under the 
isomorphism ^ has the property 

^(0)(x)i = 1 V X G 
we call this monomial function the one function 1 := ^'(0). 

Definition 17 (Notational Definition) A directed graph 

G = {Vg, Eg, ttg-.Eg^Vgx Vg) 

that allows self loops and parallel directed edges is called digraph. 

Definition 18 Let M be a nonempty finite set. Furthermore, let n := \M\ be the cardinality of M. 
A numeration of the elements of M is a bijective mapping 

/:M^{l,..,n} 

Given a numeration f of the set M we write 

M = {/!,...,/„} 

where the unique element x G M with the property f[x) = i G {1, is denoted as fi. 

Definition 19 (Notational Definition) Let f G MF^(Fq) be a monomial dynamical system 
and G = {Vg, Eg, ttg) digraph with vertex set Vg of cardinality \Vg\ = n. Furthermore, let 
F := "if^^i^f) be the corresponding matrix of f. The digraph G is called dependency graph of f iff 
a numeration a: M ^ {!,..., n} of the elements ofVG exists such that\/ i,j G {l,...,n} there are 
exactly Fij directed edges Oj aj in the set Eg, i.e. 



TTf''{{ai,aj)) 



Fij 



Remark 20 It is easy to show that if G and H are dependency graphs of f then G and H are 
isomorphic. In this sense we speak from the dependency graph of f and denote it by Gf = {Vf, Ef, 
TTf). Our definition of dependency graph differs slightly from the definition used in fd^. 

Definition 21 (Notational Definition) Let G = {Vg, Eg, ttg) be a digraph. Two vertices a,b ^ 
Vg are called connected if there is a t G Nq and (not necessarily different) vertices vi, ...,Vt G Vg 
such that 

a ^ Vl ^ V2 ^ ... ^ Vt —> b 

In this situation we write a -w^ b, where s is the number of directed edges involved in the sequence 
from a to b (in this case s = t+1). Two sequences a b of the same length are considered different 
if the directed edges involved are different or the order at which they appear is different, even if the 
visited vertices are the same. As a convention, a single vertex a ^ Vg is always connected to itself 
a o by an empty sequence of length 0. 

Definition 22 (Notational Definition) Let G = {Vg, Eg, t^g) be a digraph and a,b Vg two 

vertices. A sequence a b 

a vi ^ V2 ... ^ vt b 

is called a path, if no vertex Vi is visited more than once. If a = b, but no other vertex is visited 
more than once, a b is called a closed path. 
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Definition 23 Let G = (Vg, Eq, ttg) be a digraph. Two vertices a,b G Vq are called strongly 
connected if there are natural numbers s, t € N such that 



a b and b -^t a 



In this situation we write a ^ b. 



Theorem 24 (and Definition) Let G = {Vg, Eq, t^g) be a digraph. ^ is an equivalence relation 
on Vg called strong equivalence. The equivalence class of any vertex a G Vg is called a strongly 
connected component and denoted by *li C Vg- 

Proof. See Definition 3.1 (2) in [6]. ■ 

Definition 25 Let G = (Vg) Eg, t^g) be a digraph and a € Vg one of its vertices. The strongly 
connected component *7i C Vg is called trivial iff*li= {a} and there is no edge a ^ a in Eg- 

Definition 26 Let G = {Vg, Eg, ttg) be a digraph with vertex set Vg of cardinality \Vg\ = n and 
Vg = {ai, ...,a„} a numeration of the elements ofVG. The matrix A G M{n x n; No) whose entries 
are defined as 

Aij := number of edges ai aj contained in Eg 
for i,j = 1, ...,n is called adjacency matrix of G with the numeration a. 

Theorem 27 Let G = {Vg, Eg, t^g) be a digraph with vertex set Vg of cardinality \Vg\ = n and 
Vg = {ai,...,o„} a numeration of the elements ofVG. Furthermore, let A G M{n x n; No) be its 
adjacency matrix (with the numeration aJ, m G N a natural number and 

B ■.= A"^ £ M{n X n;No) 

the mth power of A. Then y i,j G {1, n} the entry Bij of B is equal to the number of different 
sequences a, aj of length m. 

Proof. The proof of this well-known result can be found in [TO]. ■ 

Remark 28 Let f G MF^{Fq) be a monomial dynamical system. Furthermore, let Gf = (Vj, Ef, 
TTf) the dependency graph of f and Vf = {ai, ...,a„} the associated numeration of the elements of 
Vf. Then, according to the definition of dependency graph, F := ^'^^(/) (the corresponding matrix 
of f) is precisely the adjacency matrix of Gj with the numeration a. Now, by Remarks[T^ and 64 
we can conclude 

^r-i(/™) = mre(iq(F"^) (5) 



3 Characterization of fixed point systems 

The results proved in the previous section allow us to link topological properties of the dependency 
graph with the dynamics of /. We will exploit this feature in this subsection to prove some char- 
acterizations of fixed point systems stated in terms of connectedness properties of the dependency 
graph. At the end of this section we also provide a more algebraic sufficient condition. 

Theorem 29 Let Fq be a finite field and f G MF^iYq) a monomial dynamical system. Then f is 
a fixed point system with (1, 1)* G as its only fixed point if and only if its dependency graph 
only contains trivial strongly connected components . 
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Proof. By Remark I28| F := ^'^^(/) is the adjacency matrix of the dependency graph of /. If the 
dependency graph does not contain any nontrivial strongly connected components, every sequence 
a b between two arbitrary vertices can be at most of length n — 1. (A sequence that revisits 
a vertex would contain a closed sequence, which is strongly connected.) Therefore, by theorem [271 
3 m € N with m < n such that = (the zero matrix in M{n x n; Nq)). Now, according to 
equation ^ we have 

^i-^iD = mredqiF"^) = mredq{0) = 

and consequently 

^-\f) = V r > m 

Thus 

/ '" = 1 V r > m 
If, on the other hand, there is an m € N such that 

jm+X = fm 

applying the isomorphism (see Remark USD we obtain 

^■(m+A) = ^-m = V A € N 

and (see equation ([5])) 

mredg(F"+^) = mredgiF"") = V A € N 
It follows from equation (jlj) (See also Remark I64p 

pm+a = V a e No 

Now by theorem 1271 there are no sequences a b between any two arbitrary vertices a, b of length 
larger than m — 1. As a consequence, there cannot be any nontrivial strongly connected components 
in the dependency graph of /. ■ 

Definition 30 A monomial dynamical system f € MF^(Fq) whose dependency graph contains 
nontrivial strongly connected components is called coupled monomial dynamical system. 

Definition 31 Let G = (Vg, Eq, t^g) be a digraph, m € N a natural number and a, 6 € Vg two 
vertices. The number of different sequences of length m from a to b is denoted by Smia, 6) G Nq . 

Remark 32 Let G = {Vg, Eg, ttg) be a digraph with vertex set Vq of cardinality n := \Vg\ and 
Vg = {ai,...,a„} a numeration of the elements ofVG- Furthermore, letm^N be a natural number 
and A € M{n x n; Nq) the adjacency matrix of G with the numeration a. Then by Theorem 21 we 
have 

Sm{o-i, Oj) = [A )ij 

Theorem 33 Let Fq be a finite field, f G MF^(Fq) a coupled monomial dynamical system and 
Gf = (Vf, Ef, 11 f) its dependency graph. Then f is a fixed point system if and only if there is an 
m € N such that the following two conditions hold 

1. For every pair of nodes a,b Vf with a -w^ b there exists for every A € N an € Z such 
that Sm+x{a, b) = Sm{a, b) + ax{q - 1) / 0. 

2. For every pair of nodes a,b £ Vj with Sm(o, b) = it holds Sm+xia, 6) = V A G N. 
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Proof. Let Vf = {ai, .■■,an} be the numeration of the vertices. If / is a fixed point system, 3 
m G N such that 

By applying the homomorphism ^'"^ we get (see Remark I ISp 

^•(m+A) = i?-*" V A e N (6) 

By Remark 1281 it follows 

mredg(F'^+^) = mredq{F'^) V A G N 

Let i,j G {1, ...,n}. If, on the one hand, (F ™)^^ = then by © we would have = V 

A G N. Consequently, by 2. of Lemma [6] we have 

= V a G No 

Now by theorem 1271 there are no sequences Oj Uj of length larger than m — 1. In other words, 
2. follows. If, on the other hand, {F-"')ij / then by ([6]) we would have {F<"^+^^)ij = {F-'^)ij / 
V A G N. Consequently, by 2. and 4. of Lemma [6] 3 ax G Z such that 

{F^+% = {F^h, + axiq - 1) V A G N 

In other words 1. follows. To show the converse we start from the following fact: Given 1. and 2. 
and according to Theorem [27] and Remark [28] 

If {F"')ij = 0, then = (F™),, V A G N 

and 

if / 0, then 3 OA G Z : {F"'+% = {F'^)^J + ax{q -1)/0VAgN 

Now by 2. and 4. of Lemma [U] we have 

mredgiF""-^^) = mredgiF"") V A G N 

and by [28] 

Thus, after applying the isomorphism ^ 

■ 

The following parameter for digraphs was introduced by 
Definition 34 Let G = (Vg, Eg, ttg) be a digraph and a G Vg one of its vertices. The number 

C(a) := min \u — v\ 

is called the loop number of a. If there is no sequence of positive length from a to a, then C{a) is 
set to zero. 
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Lemma 35 (and Definition) Let G = (Vg, Eq, ttg) be a digraph and a € one of its vertices. 
If^is nontrivial then for every b £ ^ it holds 

£{b) = £{a) 

Therefore, we introduce the loop number of strongly connected components as 

C^Tt) := C{a) 
Proof. See the proof of Lemma 4.2 in [6]. ■ 

Remark 36 The loop number of any trivial strongly connected component is, due to the convention 
made in the definition of loop number, equal to zero. 

Corollary 37 Let Fg be a finite field, f E MF"(Fg) a coupled monomial dynamical system and 
Gf = {Vf, Ef, TTf) its dependency graph. If f is a fixed point system then the loop number of each 
of its nontrivial strongly connected components is equal to 1. 

Proof. Let m S N be as in the statement of the previous theorem. Let ^ C V/ be a nontrivial 
strongly connected component. For every 6 G ^ we have that b is strongly connected with itself. 
Therefore, for every s G N there is a t > s such that b —^t b. In particular, there must be a u € N 
with u > m such that b b, i.e. Su(b,b) > 1. By 2. of the previous theorem we know that 
Smib, b) 7^ 0, otherwise Su(6, b) = 0. Now from 1. of the previous theorem we know 

3axeZ: Sm+x{b, b) = s^{b, b) + ax{q -1)/0VA€N 

and in particular 

Sm+x{b,b)^Oy XeN 
Therefore, V A G N there are sequences b -^m+x b. Thus >C(Y) = C{b) = 1. ■ 

Definition 38 Let G = (Vg, Eq, ttg) be a digraph and a,b £ Vg two vertices. The vertex a is 
called recurrently connected to b, if for every s € N there is a u > s such that a b. 

Lemma 39 (and Definition) Let G = (Vg, Eg, ttg) be a digraph with vertex setVG of cardinality 
IT- '■= \^g\- Two vertices a,b € Vg are connected through a sequence a b of length t > n — 1 if 
and only if a is recurrently connected to b. 

Proof. If there is a sequence a b of length t > n — 1, then it necessarily revisits one of its 
vertices, in other words, there is a c G Vg such that 

a -^t 6 = a ^ ... — > c — > ... — 5- c — >■ ... — > 6 

Now a sequence a b can be constructed that repeats the loop around c as many times as desired. 
The converse follows immediately from the definition of recurrent connectedness. ■ 

Remark 40 Let G = {Vg, Eg, ttg) be a digraph with vertex set Vg of cardinality n := \Vg\. Then 
for any two vertices a,b £ Vg it holds: Either a is recurrently connected to b or there is an m ^'H 
with m < n such that no sequence a -^t b of length t > m exists. 

Lemma 41 Let G = (Vg, Eg, t^g) be a digraph and U Vg a nontrivial strongly connected 
component. Furthermore, let t := C{U) be the loop number ofU. Then for each a,b (z U there is 
an m € N such that the graph G contains sequences a ~^m+xt b of length m + At V A G N. 
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Proof. See the proof of Proposition 4.5 in [6]. ■ 

Theorem 42 Let G = {Vg, Eq, ttg) be a digraph containing nontrivial strongly connected compo- 
nents. If the loop number of every nontrivial strongly connected component is equal to 1 then there 
is an m (zf^ such that any pair of vertices Oj , aj € Vg with ai recurrently connected to aj satisfies 

Sm+\{ai,aj) > V A € No 

Proof. Let Vg = {ai, fln} be the numeration of the vertices and ai, aj £ Vg- ai is recurrently 
connected to aj, then necessarily there is a sequence o-j that visits a vertex contained in a 

nontrivial strongly connected component. In other words, 3 at € V/- and a sequence ai -^s CLj such 
that Ofc is nontrivial and 

Oj 0.j = Ci —>■ ... ttk ^ ... — > Oj 

By Lemma HH there is a G N such that there are sequences '^m^+x V A G Nq. Now V 
A E No we can construct a sequence 

Now, if we consider among all pairs i,j £ {1, ...,n} such that € Vg is recurrently connected to 
flj £ Vg the maximum m of all values we can state: 3 m € N such that any pair of recurrently 
connected vertices ai , aj G Vg satisfies 

Sm+xiai,aj) > V A G No 



Theorem 43 Let F2 be the finite field with two elements, f £ MF^ {F2) a Boolean coupled 
monomial dynamical system and Gf = (Vf, Ef, tt f) its dependency graph, f is a fixed point system 
if and only if the loop number of each nontrivial strongly connected components of Gf is equal to 1 . 

Proof. The necessity follows from Corollary [371 Now assume that each nontrivial strongly 
connected components of Gf has loop number 1 and let Vf = {ai,...,a„} be the numeration of 
the vertices. Furthermore let F := ^~^(/) be the corresponding matrix and consider vertices 
ai,aj G Vf. By Remark HOl either Oj is recurrently connected to Oj or there is an ito G N with 
uq <n such that no sequence at aj of length t > uq exists. If the latter is the case, then 

(F"o+^)ij = V A G No 

On the other hand, if Oj is recurrently connected to aj, then by Theorem 1421 there is an tuq G N 
such that 

/ V A G No 

Therefore, we have for m := max(mo,tio) that 

{F"'+% / V A G No or {F"'+% = V A G No 
Summarizing we have by 2. of Lemma [6] 

mredg(F"^+^) = mredqiF'^) V A G N 



and by [28] 

pim+X) ^ jp m V A G N 
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Thus, after applying the isomorphism ^ 

jm+X = V A G N 

■ 

Remark 44 The statements of the previous theorems together with the Remark^ about zero func- 
tions as components constitute the statement of Theorem 6.1 in J^. 

In the following two corollaries we provide alternative proofs to the claims made in Corollary 
6.3 and Theorem 6.5 of j6j: 

Corollary 45 (and Definition) Let F2 the finite field with two elements and f G MF^ (F2) the 
coupled monomial dynamical system defined by 

fi{x) = x^" 

i-l 

where aij G Eg, i = 1, ...,n, j = 1, ...,i — 1. Such a system is called a Boolean triangular system. 
Boolean triangular systems are always fixed point systems. 

Proof. From the structure of / it is easy to see that every strongly connected component of the 
dependency graph of / is either trivial or has loop number 1. ■ 

Corollary 46 Let F2 the finite field with two elements, f G MF^ (F2) a fixed point system and 
j,i G {1, ...,n}. Consider the system g G MF^(F2) defined as gk{x) := fk{x) V A; G {1, ...,n}\j and 
gj{x) := Xifj{x) V X G F2. Then g is a fixed point system if there is no sequence o.j from 

to Oj or if*^ or'oj are nontrivial. 

Proof. If i = j then Eg contains the self loop Oj — > Oj and ^ becomes nontrivial (if it wasn't 
already) with loop number 1. li i ^ j then we have two cases: If there is no sequence a, Oj, 
then adding the edge aj — > Oj (which might be already there) doesn't affect Oj 7^ aj . If there is a 
sequence Oj Oj then adding the edge aj — > Oj (which might be already there) forces o^' = o^. 
Now since by hypothesis Oj or aj are nontrivial and / is a fixed point system, then 

£( V) = C-i'oj) = 1 



Definition 47 Let Fg be a finite field, f G MF^(Fq) a monomial dynamical system and Gf = {Vj, 
Ef, TTf) its dependency graph, f is called a (q — \)-fold redundant monomial system if there is an 
G N such that for any pair a,b ^ Vf with a recurrently connected to b, the following holds: 

y m> N 3 Oabm e No : s.m(a, b) = aabm{q - 1) 

Lemma 48 Let Fg be a finite field, f G MF^(Fq) a coupled (q — \)-fold redundant monomial 
dynamical system and Gf = (Vj, Ej, ttj) its dependency graph. Then f is a fixed point system if 
the loop number of each nontrivial strongly connected component of Gf is equal to 1. 
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Proof. Let Vf = {ai,...,a„} be the numeration of the vertices and F := ^'~^(/) be the cor- 
responding matrix of /. Consider two arbitrary vertices ai,aj € Vf. By Remark 140 ( either Oj is 
recurrently connected to Uj or there is an m-o € N with uiq < n such that no sequence a b of 
length t > niQ exists. If the latter is the case, then 

^Fm>+\^.. = V A G No 

On the other hand, if a, is recurrently connected to aj, then by Theorem 1421 there is an mi € N 
such that 

(a„a,) >0V7G No (7) 

Consider now m2 := max(n,mi). Due to the universality of mi in the expression ([7]), for any pair 
of vertices ai,aj € Vq with recurrently connected to aj there is a sequence Oj -^^^2+7 o.j of length 
m2 + 7, in particular S(m2+7)(^j) "^i) > ^ 7 ^ I^o- Now, let N be the constant in Definition 1471 and 
ms := max(A^, 7712). Now, by hypothesis, 3 Oij-y G N such that 

5(^3+7) (aj,aj) = aij^{q - 1) V 7 G No 

Thus 

s^rn3+-y)iai,aj) = aij^{q - 1) = aijo{q - 1) + {oij^ - aijo)iq - 1) 
= Sm3{ai,aj) + (oij^ - aijo){q - 1) V 7 G No 

Summarizing, since mo < n < m2 < m^, we can say V G {1, ■■■,n}, depending on whether Oj 
and aj are recurrently connected or not, 

(^"^3+^) = V A G No 

or 

3axeZ: {F'^-'+% = {F^^ij + ax{q - 1) V A G No 
Now, by 2. and 4. of Lemma [6] it follows 

mredg(F'"3+^) = mredg{F"^^) V A G N 

and by [28] 

pim3+\) ^ p ma V A G N 

Thus, after applying the isomorphism ^ 

jms+A = /"^S V A G N 



Theorem 49 Let Fg be a finite field, f G MF"(Fq) a coupled monomial dynamical system and 
Gf = (Vf, Ef, TTf) its dependency graph. Then f is a fixed point system if the following properties 
hold 

1. The loop number of each nontrivial strongly connected component of Gf is equal to 1. 

2. For each nontrivial strongly connected component ^ '^Vf and arbitrary b,cG ^ , 

si(6,c) / 0^si(6,c) =q-l 
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Proof. Let Vj = {oi, ...,an} be the numeration of the vertices and F := ^'^^(/) be the corre- 
sponding matrix of /. Consider two vertices ai,aj G Vf such that is recurrently connected to aj. 
Then by Theorem there is an rn-i G N such that 

Smi+'yiai, aj) > V 7 G Nq (8) 

Consider now m2 := max(n,mi). Due to the universality of mi in the expression ([8]), for any pair 
of vertices aj, aj G Vq with recurrently connected to aj there is a sequence Oj -^^^2+7 ^^j of length 
1712 + 7. Since m2 + 7 > n — 1, necessarily 3 ak^,ai^ G Sfc^ such that a^^ is nontrivial and 

«i "^(012+7) '^i = o-i—" ■■■ «fc-y ""^t % ■■■ — *• «i (9) 

{t depends on i,j and 7). Now, by hypothesis, every two directly connected vertices a,b & oj^ 
are directly connected by exactly q — 1 directed edges. Therefore, for any sequence ak^ -^t ai^ 
of length t G N there are {q — 1)* different copies of it and we can conclude3 a G N such that 
stio-kj , ai_^) = a{q — 1). As a consequence, there are a{q — 1) different copies of the sequence ([9]). 
Since we are dealing with an arbitrary sequence "^(012+7) "^j of fixed length m2 + 7, 7 G No we 
can conclude that 3 Uij^ G N such that 

S(m2+^)(ai,aj) = aij^{q - 1) V 7 G No 

Thus / is a coupled {q — l)-fold redundant monomial dynamical system and the claim follows from 
Lemma HSJ ■ 

Corollary 50 Let F2 he the finite field with two elements, f G MF^ (F2) a Boolean monomial 
dynamical system and F := ^'^^(/) G M(n x n; E2) its corresponding matrix. Furthermore, let 
Fq be a finite field and g G MF^(Fq) the monomial dynamical system whose corresponding matrix 
G := ^~^{g) G M(n x n; Eg) satisfies^ i,j G {l,...,n} 

_(q-lifFij = l 

If f is a fixed point system then g is a fixed point system too. 

Proof. Let Gf = {Vf, Ej, vrj) be the dependency graph of /. By the definition of g, one can 
easily see that the dependency graph Gg = {Vg, Eg, -Kg) of g can be generated from Gf by adding 
q — 2 identical parallel edges for every existing edge. Obviously Gf and Gg have the same strongly 
connected components. If doesn't contain any nontrivial strongly connected components, then 
Gg wouldn't contain any either and by Theorem [29] g would be a fixed point system. If, on the 
other hand, Gf does contain nontrivial strongly connected components, then by Theorem 1431 each 
of those components would have loop number 1. From the definition of g it also follows for any pair 
of vertices a,b ^ Eg 

si{a, 6) / =^ si(a, 6) = g — 1 
By the previous theorem g would be a fixed point system. ■ 

Example 51 (and Corollary) Let Fq be a finite field and f G MF^{Fq) the coupled monomial 
dynamical system defined by 

fi{x) = xl'^ 

i-l 
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where Uij G Eg, i = 1, ...,n, j = 1, — 1 are not further specified exponents. Such a system is 
called triangular. It is easy to see that the dependency graph of f contains n one vertex nontrivial 
strongly connected components. Each of them has a {q — \)- fold self loop. Therefore, hy the previous 
Theorem, f must he a fixed point system. 

Theorem 52 Let Yq he a finite field, f € MF^(Fq) a coupled monomial dynamical system and 
Gf = {Vf, Ef, TTf) its dependency graph. Then f is a fixed point system if for every vertex a € Vf 
that is recurrently connected to some other vertex b G Vf the edge a ^ a appears exactly q — 1 times 
in Ef, i.e. 



7r^^((a,a)) 



1 



Proof. Let Vf = {ai, ...,a„} be the numeration of the vertices and F := ^'^^(/) be the corre- 
sponding matrix of /. Consider two vertices ai,aj € Vf such that ai is recurrently connected to aj. 
Then by Theorem |32] there is an mi G N such that 



'mi+7 



(a„Oj-) > V7 GNo (10) 



Consider now m2 := max(n, mi). Due to the universality of mi in the expression (jlOp . for any 
pair of vertices ai,aj G Vq with Oj recurrently connected to aj there is a sequence Oj -^^2+7 
of length 1712 + 7- Consider one particular sequence ~-^m2+'y o.j of length m2 + 7 and call it 
w-y := ai ~-!-m2+7 hypothesis there are exactly q — 1 directed edges Oj — > Oj. Therefore, there 

are q — I copies of the sequence Wj. Since we are dealing with an arbitrary sequence -^(^2+7) 
of fixed length m2 + 7, 7 € No we can conclude that 3 Uij^ £ N such that 

•S(m2+7)(«j>«i) = Oiij^iq - 1) V 7 G No 

Thus / is a coupled {q — l)-fold redundant monomial dynamical system and the claim follows from 
Lemma [ 



Example 53 (and Corollary) Let Fq he a finite field and f G ME^(Fq) a monomial dynamical 
system such that the diagonal entries of its corresponding matrix F := "^"^(f) satisfy 

Fa = q - 1 y i e {1, ...,n} 

Since every vertex satisfies the requirement of the previous theorem, f must he a fixed point system. 
This result generalizes our previous result ahout triangular monomial dynamical systems. 

We now provide a more algebraic sufficient condition for a system / G MF^{Fq) to be a fixed 
point system. 

Lemma 54 Let n G N &e a natural numher and A G M{n x n; R) a real matrix. In addition, let 
A he diagonalizahle over C. Then A^ = AV m £ N if and only i/ 3 r, s G Nq such that r + s = n 
and the characteristic polynomial charpoly{A) of A can he written as 

charpoly{A) = a{\ - IfX* 

where a G M\{0}. 

Proof. The proof of this simple linear algebraic result is left to the interested reader. ■ 



A framework to study monomial dynamical systems over a finite field 



19 



Theorem 55 Let Fg be a finite field, f G MF^(Fq) a coupled monomial dynamical system and 
F := ^~"^(/) G M{n x n; Eg) its corresponding matrix. If the matrix F (viewed as a real matrix F 
G M{n X n; N) C M{n x n; R)J has the characteristic polynomial 

charpoly{F) = a(A - 1)'A* (11) 

where a G Z\{0}, r, s G No such that r + s = n and the geometric multiplicity of the eigenvalues 
and 1 is equal to the corresponding algebraic multiplicity, then f is a fixed point system. 

Proof. It is a well-known linear algebraic result that if there is a basis of eigenvectors of a matrix, 
the matrix is diagonalizable. By the hypothesis this is the case for F. Therefore, by the previous 
Lemma 

= F V m G N 

Now, by Remarks [15] and [Ml we consequently have V m G N 

^r-i(/™) = = mredgiF"") = mredq{F) = F 
After applying the isomorphism ^ we get 

/™ = / VmGN 



Remark 56 Let Fg be a finite field, f G MF'^iFg) a coupled monomial dynamical system and 
F := ^~^{f) G M(n x n; Eg) its corresponding matrix. The matrix F viewed as the adjacency 
matrix of the dependency graph Gj = (Vf, Ej, ttj) of f satisfies 

F"" = F V m G N 

if and only if for each pair of vertices a, 6 G V/ the value Smio-,b) is constant for all m G N. In 
other words, a and b are either disconnected or for every length m G N they are connected with the 
same degree of redundancy. 

Example 57 Consider the monomial system g G MF^ (F3) defined by the matrix 

/I 1 0\ 

110 

G:= 10 

1 

\0 0/ 



It is easy to show that 

charpoly{G) = (A - 1)^X'^ 
However, g is not a fixed point system. This shows that the condition 177]) alone is not .sufficient. 
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4 An algorithm of polynomial complexity to identify fixed point systems 
4.1 Some basic considerations 

Definition 58 Let X he a nonempty finite set, n G N a natural number and f : X^ a time 

discrete finite dynamical system. The phase space of f is the digraph with node set X", arrow set 
E defined as 

i?:={(x,y) GX" xX- | f{x)=y} 

and vertex mapping 

vr : ^ ^ X" X X" 

(x,y) ^ {x,y) 

Remark 59 Due to the finiteness of X it is obvious that the trajectory 

X, fix), f^{x),... 

of any point x S X" contains at most jX"! = |X|" different points and therefore becomes either 
cyclic or converges to a single point y € X with the property f{y) = y (i.e. a fixed point of f). 
Thus, the phase space consists of closed paths of different lengths between 1 (i.e. loops centered on 
fixed points) and \X^\ = |X|" and directed trees that end each one at exactly one closed path. The 
nodes in the directed trees correspond to transient states of the system. 

According to our definition of monomial dynamical system / G MF^(Fq), the possibility that 
one of the functions fi is equal to the zero function is excluded (see Definition [3] and Remark 
H]). Therefore, the following algorithm is designed for such systems. However, in this algorithmic 
framework it would be convenient to include the more general case (as defined in [5J and |6j), i.e. 
the case when some of the functions fi can indeed be equal to the zero function. In the vein of 
Remark U this actually only requires some type of preprocessing. The preprocessing algorithm will 
be described and analyzed in the Appendix. 

Our algorithm is based on the following observation made by Dr. Michael Shapiro about 
general time discrete finite dynamical systems: By the previous remark, a chain of transient states 
in the phase space of a time discrete finite dynamical system / : X" X" can contain at most 
s := |X"| — 1 = |X|" — 1 transient elements. Therefore, to determine whether a system is a fixed 
point system, it is sufficient to establish whether the mappings and f^~^^ are identical for any 
r > s. In the case of a monomial system / G MF^(Fq), due to Theorem 1141 we only need to look 
at the corresponding matrices F'^, g M(n x n;Eq). Computationally it is more convenient 

to generate the following sequence of powers 

(F-2)-2 = i7-4^ (^■4y2^^.8^ [E^y^ = F^^ . . . , F<^'^ 

To achieve the " safe" number of iterations | | — 1 = — 1 we need to make sure 

2* > - 1 

This is equivalent to 

t > log2(g" - 1) 
To obtain a natural number we use the ceil function 

t := cei/(log2(g" - 1)) (12) 
Thus we have, due to the monotonicity of the log function, 

t < log2((?" - 1) + 1 < log2(<7") + 1 = nlog^iq) + 1 
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4.2 The algorithm and its complexity analysis 

The algorithm is fairly simple: Given a monomial system / € MF^(Fq) and its corresponding 
matrix F := ^"^(/) G M(n x n;Eq) 

1. With t as defined above (|12|) . calculate the matrices A := F'"^^ and B := FA. This step 
requires t + 1 matrix multiplications. 

2. Compare the entries Aij and Bij. This step requires at most comparisons. (This 
maximal value is needed in the case that / is a fixed point system). 

3. / is a fixed point system if and only if the matrices A and B are equal. 

It is well known that matrix multiplication requires 2n^ — v? addition or multiplication oper- 
ations. Since t + 1 < n\og2{q) + 2, the number of operations required in step 1 is bounded above 
by 

(2n^ -n2)(nlog2(g) + 2) 

Summarizing, we have the following upper bound N{n, q) for the number of operations in steps 1 
and 2 

N{n,q) := (2n^ - n^)(n log2(g) + 2) + 
For a fixed size q of the finite field Fq used it holds 

hm ^ = 21og2(g) 

and we can conclude N{n, q) G 0{ n^) for a fixed q. The asymptotic behavior for a growing number 
of variables and growing number of field elements is described by 

hm ^- = 2 

Thus, N{n,q) G 0{ n^log2(g)) for n, g — > oo. 

It is pertinent to comment on the arithmetic operations performed during the matrix multiplica- 
tions. Since the matrices are elements of the matrix monoid M(n x n; Eg), the arithmetic operations 
are operations in the monoid Eg. By the Lemmas I61|and 1601 the addition resp. the multiplication 
operation on Eg requires an integer number additioro resp. multiplication and a reduction as de- 
fined in Lemma El The reduction redg{a) of an integer number a G Nq, a > g is obtained as the 
degree of the remainder of the polynomial division r" ( r'^ — r). According to 4.6.5 of [13] this 
division requires 

0(2(deg(r«) - deg(r« - r))) = 0(2(a - q)) 

integer number operations. However, we know that the reductions redg{.) are applied to the result 
of (regular integer) addition or multiplication of elements of Eg and therefore 

2{q-l)-q = q-2 
(q-l)2 _q = q2 _q^l 

As a consequence, in the worst case scenario, one addition resp. multiplication in the monoid Eg 
requires 0{q) resp. 0{q^) regular integer number operations. 

Since Eg is a finite set and only the results of pairwise additions and pairwise multiplications 
are needed, while the algorithm is running, these numbers are of course stored in a table after the 
first time they are calculated. 



< 



^ See Chapter 4 of [IJ' for a detailed description of integer number representation and arithmetic in typical computer 
algebra systems. 
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A Appendix 

A.l Some simple algebraic results 

Lemma 60 Let Fg be a finite field and a, 5 E No nonnegative integers. Then it holds 

redq{ab) = redq{redq{a)redq{h)) 

Proof. We have V x e 
and by Lemma [6] 

redq{ab) = redq{redq{a)redq{h)) 

■ 

Lemma 61 Let Yq he a finite field and a, 6 G No nonnegative integers. Then it holds 

redq{a + h) = redq{redq{a) + redq{b)) 

Proof. By the division algorithm 3i Qa , gt , Qa+b , ra , ri, , ra+b G Fg[r] such that 

r'^ = gair'^ - r) + ra = gair'^ - r) + T-^'""'^-^ 
r' = gb{r'^-r)+rk = gb{T''-T)+T''''^^''^ 

From the first two equations follows 

T^^" = gagbir" - rf + garb{T^ - r) + rMr"^ - r) + r'-dq(a)+redq{V) 
Applying the division algorithm to 'T-'r'^dq{a)+Tedq{b) ^ ^ Fg[T] such that 

t'^+' = gagb{r''-Tf+garb{T''-T)+ragb{r''-r)+gr{T''-r)+rr 

= {gagbir'^ - r) + gar^ + r^gb + 5r)(r'? - r) + r'--Mredq{a)+redqib)) 

From the uniqueness of quotient and remainder it follows 

^redq(a+b) _ ^redq{redq{a)+redq{b)) 

and consequently 

redq{a + b) = redq{redq{a) + redq{b)) 
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Theorem 62 (and Definition) Let Fg be a finite field. The set 

Eg = {0,l,...,{q-2),{q-l)} CZ 

together with the operations of addition a(Bb := redq{a + b) and multiplication a»b := redq{ab) is a 
commutative semiring with identity 1. We call this commutative semiring the exponents semiring 
of the field Fg. 

Proof. First we show that Eg is a commutative monoid with respect to the addition ©. The 
reduction modulo the ideal (r'^ — r) ensures that Eg is closed under this operation. Additive 
commutativity follows trivially from the definition. The associativity is easily shown using Lemma 
inland the fact that c € -Eg c = redg{c). It is trivial to see that is the additive identity element. 
Eq is also a commutative monoid with respect to the multiplication • : The reduction modulo the 
ideal {t'^ — r) ensures that Eg is closed under this operation. Multiplicative commutativity as well 
as the fact that 1 is the multiplicative identity follow trivially from the definition. The associativity 
is shown using Lemma [60] and the fact that c £ Eq <^ c = redq{c). The proof of the distributivity 
is a straightforward verification. ■ 

Lemma 63 Let n G N 6e a natural number, Fg be a finite field and Eq the exponents semiring of 
Fg. The set M{n x n; -Eg) ofnxn quadratic matrices with entries in the semiring Eg together with 
the operation ■ of matrix multiplication over Eg is a monoid. 

Proof. The matrix multiplication ■ is defined in terms of the operations © and • on the matrix 
entries, therefore M{n x n; Eg) is closed under multiplication. The proof of the associativity is a 
tedious but straightforward verification. The identity element is obviously the unit matrix /. ■ 

Remark 64 Since the entries for the matrix product D = A ■ B are defined as 

Dij = Ail • Bij © Ai2 • © ... © Ain • Bnj 

according to the definitions of the operations • and © we can write 

Dij = redqiAiiBij) © redg{Ai2B2j) © ... © redg{AinBnj) 

= redq{redq{AiiBij) + redq{Ai2B2j) + ... + redq{AinBnj)) 

Now, by Lemma\6^ we have 

Dij = redg{AiiBij + Ai2B2j + ... + AinBnj) 

As a consequence, if we define the following reduction operation for matrices with nonnegative 
integer entries 

mredg : M{n x n; Nq) M{n x n; Eg) 
Aij I > redg{Aij) 

then the following property holds for U,V M{n x n; No) and W := UV € M{n x n; No) 

mredgiUV) = mredg(U) ■ mredq{V) 

It can be easily shown that M{n x n; Nq) is a monoid and mredq : M{n x n; Nq) M{n x n;Eq) 
a monoid homomorphism. In addition, by 2. of Lemma we can conclude 



mredq{A) = 0^ A = 



(13) 
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A. 2 The preprocessing algorithm 

We start with the definition of monomial dynamical system according to [6] and [S]: 

Definition 65 Let Yq he a finite field. A map f : — > is called a monomial dynamical 
system over Fq if for every i G {l,...,n} there exists a tuple {Fn, Fin) G Eq and an element 
Oi € {0, 1} C Fq such that 

= aixf V X G F^ 

In order to use the algorithm described in Section 4 to determine whether such a monomial 
dynamical system is a fixed point system we need to preprocess the system in the sense of Remark 
m To accomplish this task algorithmically, we add an element —oo to our exponents semiring Eq^ 
(See Definition [2] and Theorem 1621): 

E~q := Eq U{-oo} 

The arithmetic with this new element is as follows 

a © — oo = — oo © a = — cxD V a € -Eg 
a • —oo = —oo • a = — oo V a S £'g\{0} 
O.a = a»0 = OVaG^ 

The addition is due to the additive "absorption property" of — oo obviously associative. The same 
holds for the multiplication, since both and — oo show the multiplicative "absorption property" 
(although wins over — oo). With this rules we are already able to multiply pairs of matrices with 
entries in Eq. With this extended exponents set we can represent the monomial dynamical systems 
defined above as follows: 

Definition 66 Let Fq be a finite field. A map f : F^ — > F^ is called a monomial dynamical system 
over Fq if for every i G {f,...,n} there exists a tuple {En, Fin) £ Eq or a tupel {En, Fin) G 
{— oo}" such that 

fi{x)=X^'\..X^^- \f XeF^q 

Now we describe the preprocessing algorithm: Given a monomial system f G ME^(Fq) and its 
representing matrix F G M{n x n; Eg) 

1. Initialize Li := and L2 := and an array v of length n to zero. 

2. For k from I to n do L2 := L2 + 1 and v[k] := 1 if and only if E^i = —00. 

3. Compare Li and L2. If Li = L2 or L2 = n, construct the matrix 

F' G M((n-L2) X (n-Ls);^^) 

by deleting the kth row and the kth column of F for all k s.t. v[k] = 1. Then return F' and 
stop. If Li < L2 and L2 < n, calculate the product F'^, set F := F'"^ as well as Li := L2 and 
go to step 2. 

4. If the returned matrix F' is the empty matrix (L2 = n) we can conclude that the system / 
is a fixed point system with (0, ...,0)* G F^ as its unique fixed point (see Remark[3]). If F' is 
not the empty matrix, the corresponding lower dimensional system /' := ^{F') needs to be 
analyzed with the algorithm described in Section 4. 
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Step 1 implies n + 2 initializations. Step 2 of the algorithm requires n comparisons, at most 
n additions and at most n assignments. There are 2 comparisons in step 3. Each matrix mul- 
tiplication in step 3 takes 2n^ — n? addition or multiplication operationqj] in Eq. There is one 
initialization after each matrix multiplication. The worst case scenario is given when every time 
the algorithm performs step 3, the set Li grows by one element, forcing the algorithm to perform 
n — 1 matrix multiplications. The construction of the matrix F' requires a number of comparisons 
and assignments that is obviously bounded above by 2v? . Summarizing, the worst case complexity 
of the algorithm is bounded above by 

B{n) : = (n + 2) +n(3n) +n2 + (n - l)(2n3 - + 1) + 2n^ 
= 2n^ - 3n^ + Gn^ + 4n + 1 

Since 

B(n) 
hm — = 2 

n—>oo n 

we can conclude B{n) S 0( n^). 

It is pertinent to emphasize that this preprocessing algorithm represents a primitive first attempt. 
Since the matrix multiplications dominate the complexity of the algorithm, it seems meaningful to 
try to reduce the complexity of the multiplication. Indeed, the rows with entries — oo are preserved 
during the multiplication, i.e. those rows do not need to be calculated. In addition, if the first 
element of a row in the product matrix is equal to — oo we know that all the remaining elements of 
that row are going to be equal to — oo as well. As we can see, there are possibilities of improvement. 
However, for the purposes of this paper, we are satisfied with a first working algorithm of polynomial 
complexity. 
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